Thermal conductivity of anisotropic and frustrated spin-1/2 chains 
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We analyze the thermal conductivity of anisotropic and frustrated spin-1/2 chains using analytical 
and numerical techniques. This includes mean-field theory based on the Jordan- Wigner transfor- 
mation, bosonization, and exact diagonalization of systems with A'^ < 18 sites. We present results 
for the temperature dependence of the zero-frequency weight of the conductivity for several values 
of the anisotropy A. In the gapless regime, we show that the mean-field theory compares well to 
known results and that the low-temperature limit is correctly described by bosonization. In the 
antiferromagnetic and ferromagnetic gapped regime, we analyze the temperature dependence of the 
thermal conductivity numerically. The convergence of the finite-size data is remarkably good in the 
ferromagnetic case. Finally, we apply our numerical method and mean-field theory to the frustrated 
chain where we find a good agreement of these two approaches on finite systems. Our numerical 
data do not yield evidence for a diverging thermal conductivity in the thermodynamic limit in case 
of the antiferromagnetic gapped regime of the frustrated chain. 



Introduction - Transport properties of low- 
dimensional spin systems have attracted recently inter- 
est both from the experimental and theoretical side. A 
particular motivation comes from the observation that 
magnetic excitations of one-dimensional spin systems sig- 
nificantly contribute to the thermal conductivity which 
is manifest in many exneriments on materials such as 
the spin-ladder systemtlEld (Sr,La,Ca).^4Cu2404i and the 
spin-chain compounds SrCu02 and Sr2Cu03u. Assum- 
ing elementary excitations to carry the thermal current 
and using a relaxation time ansatz for their kinetic equa- 
tion one finds extremely large mean- free paths being, for 
example, of the order of lOOOA in La5Ca9Cu2404iB. Al- 
though the magnitude of the mean-free path is currently 
an issue of intense discussion, the question arises whether 
heat transport in low-dimensional spin systems is bal- 
listic, i.e., whether intrinsic scattering of magnetic ex- 
citations is ineffective to render the thermal conductiv- 
ity finite. From the theoretical point of view this issue 
is related to the value of the so-called (thermal) Drude 
weightd -Dth which is the zero-frequency weight of the 
thermal conductivity At. A nonzero value of Z^th corre- 
sponds to a diverging thermal conductivity. This sce- 
nario is trivially realized if the energy-current operator 
is a conserved quantiiK, which is the case for the spin- 
1/2 Heisenberg chairflu. For a number of other models 
like the frustrated chain, the dimerized chains or the spin 
ladder the energy-current operator is not conserved and 
the question of nonzero Z^tii is a challenging topic. 

In this paper, we establish various numerical and an- 
alytical techniques to analyze the thermal Drude weight 
and to compute the temperature dependence of Dth (T) • 
We study the model Hamiltonian H = hi with the 
local energy-density given by 

hi = J{(5+5,- 1 + H.c.)/2 + A5f + aSi ■ S1+2}. (1) 

The XXZ model {a — 0) is integrable whereas for 



nonzero frustration the model becomes nonintegrable. 
Recently, Kliimper and SakaiQ obtained D^{T) for a = 
and < A < 1 by using the Bethe ansatzis which allows 
us to test our approaches in this regime. Exact diagonal- 
ization of finite systems uq to = 14 sites has been ap- 
plied by Alvarez and Groaa to investigate Dt\i{T) for the 
isotropic Heisenberg chain (i.e., A = 1), the frustrated 
chain, and the spin ladder. Our numerical analysis goes 
beyond this by allowing A 7^ 1 and extension to larger 
systems with A^ < 18 sites. 

Thermal conductivity - The thermal conductivity 
K is defined-by (j) = — kVT and is given by the following 
expressionll3: 



KiLo)^P dte-'^' dT{ji~t-lT)j). (2) 

Jo Jo 

j is the energy-current operatorS and f3 — 1/T is the 
inverse temperature. The current operator satisfies the 
equation of continuity: dthi — i[H,hi] — —{ji+i—ji)- For 
exchange interaction of arbitrary range, i.e., [hi±m, hi] ^ 
for m < Too, this implies 



jl = i ^ [hl-m-l,hi+n]. 
m,n— 



(3) 



In our case we have toq = 2 (see Eq. (|])) leading to 

j = "^ji = iY^[hi^2 + hi^i,hi + hi+i]. (4) 

Note that the current operator derived from Eq. (^ in- 
cludes the proper limiting form for a ^ 1 where one 
recovers the current operators of two decoupled chains. 

The conductivity k(w) may be decomposed according 
to HeK(^uj) = Dth{T) S(^uj) + Kicg(w) into a singular part 
at zero frequency and a regular part Krcg(i^)- The quan- 
tity of interest is the thermal Drude weight Dth{T) which 
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— Bethe-ansatz, Klumper et al. (2002) 

— Mean-field theory 
Low-temperature limit 



A=0,0.5,l;a=0 




FIG. 1: -Dth(r) for diflerent values of the anisotropy A = 
1, 0.5, (top to bottom) and zero frustration. Dashed hues de- 
note results obtained by Jordan- Wigner transformation and 
mean-field treatment of the interaction term. Bethe-ansatz 
results from Ref. ^ are included in the figure (thick solid 
lines). Thin dashed lines show the exact result for the low- 
temperature limit (see Eq. (bl)). 



can be computed via (see, e.g., Ref. ||) 



ZN 



\{m\3\n)[^ 



(5) 



Z is the partition function and N the number of lattice 
sites. Note that we exclusively use periodic boundary 
conditions. If j is a conserved quantity, the conductivity 
reduces to k(w) = Dt-\^{T)5{ijj) and expression simpli- 
fies to Ah = Trl3^{f )/N. 

Mean-field theory - Using a Jordan- Wigner 
transformationE3 the spin operators are mapped to spin- 
less fermionic operators cl^"*. In the case of A = and 
zero frustration a the corresponding Hamiltonian is diag- 
onal in momentum space and reads H = J^k^kC^Ck with 
a tight-binding dispersion = — Jcos(fc). A nonzero 
value of A or a leads to a four-fermion interaction term 
that can be treated approximately by Hartree-Fock (for 
details, see, e.g., Ref. ^2|), resulting in a renormaliza- 
tion of efc to efc = - J(H- 2A{A - 2a)) cos(fc). The pa- 
rameter A — ^ Jo dk cos{k)f{ik) has to be determined 
self-consistently where /(e) — l/(exp(/3e) -I- 1) is the 
Fermi function. Using Eq. (|5[) the thermal conductivity 
k{lu) — Dth{T)6{uj) can be computed directUi, Here we 
focus on the case of < A < 1 and a < acriJij, which is 
known to exhibit gapless spinonlike excitations. Results 
for Dth{T) of the XX Z model are shown in Fig. |l| (thick 
dashed lines). For comparison Bethe-ansatz resultsQ are 
included in the figure (sohd lines). In the case of A = 0, 
no approximations are necessary in the Jordan- Wigner 
approach and consequently, Jordan- Wigner and Bethe- 
ansatz results are identical. For A > 0, the main ob- 
servation is that the mean-field theory produces qualita- 
tively the right picture of the temperature dependence of 
Dth{T). Both the slope of Ah T at low temperatures 



and the position of the maximum are well predicted. De- 
viations at high temperatures are due to the neglect of 
many-particle excitations in the mean-field approxima- 
tion. 

Bosonization - In the continuum limit the physics of 
the anisotropic spin- 1/2 chain at lowjenergies is described 
by the Luttinger-liquid Hamiltonianllil 



(6) 



where </> is a bosonic field in 1 -I- 1 dimensions and Q is 
the dual field d^Q = j^drcj). K is the Luttinger parame- 
ter and V — (J7r/2)^i^ is the spinon velocity where the 
anisotropy is parametrized via A = cos (7) here. The lo- 
cal current operator j(x) is again given by the equation 
of continuity: dxj{x) — —dth{x). We obtain 



j = / dxdx(t>{x)dxQ{x). 



(7) 



The Drude weight follows from Dth — tt/?^ 
Thus we have to evaluate the two-point function 
(j(a^, r)j(0, 0)), T being the time variable. The compu- 
tation is similar to the procedure for the susceptibility 
in Ref. We change to coordinates z — vt + ix and 
z — VT — ix. By decomposing (f)(z, z) = ip{z) + (p{z) into 
its chiral parts and using the respective two-point func- 
tions such as {(p{z)(p{w)) = —{K/4:TT)\n{z — w) we obtain 



(j(x,r)j(0,0)) 



(47r)2 



1 



1 



(8) 



Before performing the space integration the imaginary 
time direction is compactified by mapping the plane (z) 
into the strip using z{() — exp(27r^//3) leading to 
the replacement vt ± ix {v[3 /tt) sin(7r """^J^" ) in Eq. 
(^). After the change of variables u = tan(7rT//9); i« — 
— i tan(i7rx/(w/3)) we finally find 



D,^{T)^—vT. 



(9) 



This coincides with Kliimper's and Sakai's analytic 
expressionQ for the low-temperature limit of the XX Z 
model if the velocity v is equal to v — ( J7r/2)2i^. How- 
ever, the result is more generally valid for models with 
the continuum limit given by the Luttinger-liquid Hamil- 
tonian. 

Exact diagonalization (ED) - In this part we 
present our results for Dti^iT) obtained by exact diago- 
nalization for finite systems with TV < 18. We start with 
the discussion of different values of the anisotropy A at 
zero frustration. Figure || shows Ah for the isotropic case 
A = 1, for a gapped, antiferromagnetic system (A — 10) 
and in the ferromagnetic regime (A = —1, —2). While we 
show in Fig. ^a) that we reproduce the results by Alvarez 
and GrosLI for system sizes of < 14, our analysis ex- 
tends this case to A^ < 18. This is due to exploiting both 
conservation of total and momentum k in the exact 
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FIG. 2: Exact diagonalization (ED) for the XXZ chain: 
(a) Isotropic chain (A = 1). Dashed(dotted) lines denote 
even(odd)-numbered systems rfor Ai' < 18 sites. Bethe-ansatz 
results by Kliimper and Sakaiu are included in the plot (thick 
solidpline). Thin solid lines show ED results by Alvarez and 
Groffl for N = 8, 10, 12, 14. (b) Antiferromagnetic, gapped 
regime (A = 10). ED for N = 8, . . . , 16 (dashed lines), 
= 18 (solid line), (c) Ferromagnetic regime (A = —1, —2). 
Note: the thermodynamic limit is reached for ~ 17. The 
insets of (b) and (c) display the exponential suppression of 
Dth(r) at low temperatures for A = 10 and A = —2 (in the 
insets, vertical axes are scaled logarithmically, horizontal axes 
reciprocally) . 



diagonalization. By comparing with the curve obtained 
from the Bethe ansatzQ (solid line in Fig. |^(a)) it can be 
seen that for a system of iV = 18 sites the thermody- 
namic Umit is reached for temperatures around T > 0.3 J 
for A = I. At low temperatures Dth{T) is exponentially 
suppressed due to the finite-size gap in the case of an 
even number of sites and divergent for an odd number. 



The latter is due to the degeneracy of the ground state 
in case of odd-numbered systems. 

For the gapped, antiferromagnetic case we choose A — 
10, shown in Fig. || (b), where the finite-size effects of 
the two-spinon gap are small. The data for Dth are con- 
vergent for T ^ 3J, but substantial finite-size effects 
are still present in the vicinity of the maximum, i.e., 
at small temperatures compared to the two-spinon gap 
8.055126jt3. At low temperatures the thermal Drude 
weight Dth (T) is expected to be exponentially suppressed 
in the thermodynamic limit. In the inset of Fig. ^ (b) 
Dth{T) is plotted logarithmically versus 1/T. If one 
fits Dth(2^) ^v-exp{—S/T) to the numerical data at low 
temperaturesEJ one finds S — 8.056J for iV = 18 sites 
with similar values of b found fcp,other N . This com- 
pares well to the two-spinon gapt3. Hence we conclude 
that mainly the elementary excitations contribute to the 
thermal conductivity at low temperatures. 

In the ferromagnetic regime (A < —1) (results are 
shown for A = — 1 and A = —2 in Fig. ||(c)) our main ob- 
servation is that convergence with iV is very good at all 
temperatures. For example, we find for A = —2 that 
the relative difference between the finite-size data for 
= 16 sites and A^ = 17 is negligibly small, namely 
\D^C^\T) - Dg-^^T)\/D^^=^^T) < 0.008 for T > 
0.05J. If one extracts (S-from a fit of r'th(2^) ~ exp(— (5/r) 
to the numerical datallZl for A — —2, we find S « 0.97J 
which coincides with the one-triplet gap —(A -I- 1) J that 
can easily be obtained from a spin-wave computation. 
The fast convergence is even more remarkable for the 
case A = — 1 where the numerical data are consistent 
with Dth(T) ~ T at low temperatures. In addition, 
there is no qualitative difference between even- and odd- 
numbered systems due to the ferromagnetic nature of the 
interaction for A < —1. 

Frustrated chain - Now we turn to the case of 
nonzero frustration. Since [H, j] ^ here, care has to 
be taken about off-diagonal matrix-elements of j if de- 
generacies occur. However, since we use classification by 
momentum k and degeneracies are lifted and do not 
play a crucial role. In Fig. ^ we show Dth{T) obtained 
numerically from Eqs. (|^) and (||) for even system sizes 
with A^ < 18 and a = 0.35, A = 1. A central result of 
this paper is that, while we observe a finite Drude weight 
at temperatures T > and all system sizes investigated, 
we still find a substantial reduction of the Drude weight 
with increasing system size at high T. This is in sharp 
contrast to the XXZ model, where finite-size effects are 
small at high temperatures (see e.g. Fig. || (a)). These 
observations clearly point to a vanishing of the Drude 
weight in the thermodynamic limit for a = 0.35. How- 
ever, the question of dissipationless thermal transport at 
arbitrary a > remains to be studied in more detail. 

Finally, we compare our mean-field approach with nu- 
merical results on finite systems and with nonzero frus- 
tration. In Fig. ^ we present the thermal Drude weight 
of systems with iV = 16 sites for a — 0,0.05,0.15 and 
A = 1, i.e., in the gapless regime. As is obvious from this 
figure, there is a good agreement between ED and the 
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FIG. 3: Thermal Drude weight Ah(r) for the frustrated 
chain with A = 1, a = 0.35 for = 8, .... 18 sites. 
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FIG. 4: Comparison between exact diagonahzation (sohd 
hnes) and mean-field approximation (dashed lines) on fi- 
nite systems with A'^ = 16 sites for various values of a = 
0, 0.05, 0.15 (top to bottom) and A = 1. 



mean-field approach regarding the temperature depen- 
dence of the thermal Drude weight. The general features 



of _Dth (r) are a reduction of the absolute value of Dth on 
increasing a, a shift of the position of the maximum to 
lower temperatures and thus a crossing of the curves for 
different a at low temperatures which are present in both 
the ED and mean- field results. Deviations at high tem- 
peratures are again understandable due to the neglect 
of many-particle excitations in the effective one-particle 
picture. 

Conclusion - We performed a detailed analysis of 
the thermal Drude weight for anisotropic and frustrated 
spin-l/2-Heisenberg chains by using mean-field theory, 
bosonization and ED. In the case of the XXZ model 
we demonstrated the applicability of these techniques 
for computing the temperature dependence of the ther- 
mal Drude weight. Using ED we obtained results on 
finite systems of < 18 sites for arbitrary values of the 
anisotropy A. In the ferromagnetic regime (A < — 1) 
of the XXZ chain the numerical data converge to the 
thermodynamic limit at arbitrary temperature for mod- 
erately small system sizes [N k. 18). Thg analytical re- 
sults compare well with the Bethe ansatzQ in the gapless 
regime of the XXZ model and to our numerics in the 
case of the frustrated chain on finite systems. Our nu- 
merical data at a = 0.35 mark a clear difference between 
the integrable XXZ case and the nonintegrable one at 
a = 0.35: while in the former case the Drude weight re- 
mains finite in the thermodynamic limit we have clear 
indications for a vanishing Drude weight at high tem- 
peratures in the latter case. Extended analysis of these 
findings for frustrated and dimerized chains and spin lad- 
ders will be the subject of a forthcoming paper. 
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